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A relativistic approach to describe nuclear and in general strongly interacting 
matter is introduced and discussed. Here, not only the nuclear forces but also 
the masses of the nucleons are generated through meson fields. Within this 
framework it is possible to calculate properties of finite nuclei at a level of 
accuracy similar to dedicated relativistic nuclear structure models. Due to the 
more general approach, a wider range of properties of hadronic states can be 
investigated. A number of results for heavy and neutron-rich nuclei toward the 
drip line are presented. 
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1. Introduction 

In modern nuclear physics one general area of research tries to push the 
boundaries of strongly interacting matter towards extremes of tempera¬ 
ture, density as well as isospin. In the case of extreme temperatures the 
experiments of choice are ultra-relativistic heavy-ion collisions that create 
a short-lived very hot fireball in the collision zone with a small net baryon 
density (i.e. about as many particles as anti-particles). At the other end 
of conditions neutron stars or more generally, compact stars, feature very 
large baryon number densities in the core and low temperatures. There¬ 
fore the analysis of neutron star observations can yield indirect hints to 
the properties of dense strongly interacting matter. However, this is not 
the only difference with respect to heavy-ion collisions. As the name indi¬ 
cates, neutron stars are predominantly made up of neutrons, if one neglects 
possible hyperonic and quark phases. Therefore the isospin density of this 
system is very high, compared to heavy-ion collisions. Here this discussion 
connects to exotic nuclei with large neutron numbers exploring the range up 
to the neutron drip line. For these systems the hadronic matter is sampled 
at around saturation density, low temperatures and relatively high isospin. 
Therefore, all these studies complement each other in the understanding of 
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strongly interacting matter under extreme conditions. 

While there are a number of theoretical approaches that try to ad¬ 
dress these various areas separately, in principle a general strong interaction 
model framework is needed that can describe all the mentioned conditions. 
This is the path we follow in this article. We describe such an approach 
for dense and hot matter that in its full version also includes quark degrees 
of freedom. However, in this article we concentrate on nuclear properties. 
The interaction parameters of the model are fine-tuned such that a good 
description of nuclei is possible. Within this framework it is then possible 
not only to study nuclei and exotic isotopes but also neutron stars and 
heavy-ion collisions opening up the possibility to establish correlations of 
observables connecting these very different regimes. 


2. Model Description 

We have followed this path of developing a general model applicable to a 
large range of densities and temperatures for a number of years, developing 
an effective chiral SU(3) hadronic (CMF) model, where the baryonic masses 
are largely generated by their coupling to scalar fields that attain a vacuum 
expectation value due to spontaneous symmetry breaking, analogously to 
simple chiral a models 1 . 

The model is described in detail in various publications 2-4 . Within 
the CMF approach baryons interact via meson exchange, where the meson 
fields are approximated by their mean field values. In consequence, the 
effective baryon masses m* read 

m* = g ia cr + g i( ;C + gitS + Snii . (1) 

where i labels the various baryons of the lowest baryonic multiplet contain¬ 
ing the nucleons as well as the A, £ + / + /°, and S 0, ~ hyperons. The scalar 
fields in this equations, a, S, correspond to the scalar quark condensates, 
a ~< uu + dd >, £ ~< ss >, and <5 ~< uu — dd >. As it is a very well es¬ 
tablished result from QCD, clearly shown in lattice gauge simulations, the 
vacuum state contains non-vanishing scalar (and isoscalar) condensates, di¬ 
rectly related to a and £. In the model approach this is achieved through 
a flavor-SU(3) self-interaction that generates non-vanishing scalar fields in 
the vacuum 4 . A mass term Srrii breaks chiral and SU(3) symmetry explic¬ 
itly. The various couplings contained in the equation result from the SU(3) 
coupling scheme. Thus, the non-zero vacuum expectation values of a and £ 
generate the baryonic masses in the vacuum, while the change of the scalar 
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fields at finite density or temperature are responsible for the scalar attrac¬ 
tion. The non-vanishing vacuum condensates are the result of the structure 
of the SU(3)-invariant scalar self-interactions, given by 3,4 

Lseif = k 0 (<J 2 + e 2 + s 2 ) + ki (a 2 + C 2 + s 2 ) 2 + k 2 ^ + 3 a 2 S 2 + 

+k 3 (a 2 - 6 2 )C + k 4 In (tT ^ )C ■ (2) 

with various coupling constants ki that have largely to be chosen to repro¬ 
duce baryon masses 4,6 . 

The vector fields ( rho° meson) and (j) are the analogous fields to 
the cr, 5, and </, respectively. The fields acquire non-zero mean fields in 
dense matter and generate repulsive forces between baryons via the linear 
interaction term 


Lbv = - ^ V’i + 9 P iP + 94>i<i>\ i’i (3) 

i 

summing over the baryons with their fields ipi- The various couplings are 
related through SU(6) relations 4 . 

Thus the interplay of the scalar and vector fields produces saturation of 
nuclear matter. In order to obtain not just reasonable infinite matter results 
but also a good description of finite nuclei, a careful fine-tuned optimisation 
of the parameters has to be performed. To illustrate the usefulness of the 
approach it is worth pointing out that the hadronic model can be readily 
extended to include quarks as well as has been discussed in'’ 8 in detail. 
Both quarks and hadrons couple to the condensate fields and therefore allow 
for smooth transition from hadronic degrees of freedom at low temperature 
and/or densities to high densities and temperatures, respectively. This can 
also be illustrated by calculating the scalar condensate a as function of 
temperature T and baryon chemical potential fib- The resulting plot is 
shown in Fig. 1. One can see the reduction of the condensate at high T 
and /if,. The model also shows a first-order liquid gas phase transition. 

On the purely hadronic side, a very recently improved fit of the coupling 
constants 6 was used to calculate all even-even nuclei up to the proton and 
neutron drip lines. The resulting errors e in binding energy compared to 
experimental data 9 amount to e = 0.17% for nuclei with baryon number 
larger than or equal to A = 50, and e = 0.12% for A < 100., which is a 
20% improvement to earlier results 4 . As it is common in the case of mean 
held approaches, larger systems can be described more accurately. 
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Fig. 1. Scalar fiel a as function of temperature and baryochemical potential. The shades 
code the value of the field. Also marked are the first-order liquid-gas phase transition as 
well as the position of the smooth deconfinement transition. 


In addition to a satisfactory description of binding energies of known nu¬ 
clei a reasonable description of surface effects are an important part of a 
successful nuclear model. In this context various deformation properties of 
nuclei were investigated. Fig. 2 shows the quadrupole deformation /?2 of 
sulfur isotopes. Also shown are measured values for /? 2 10 - The figure shows 
good agreement of data and model. Furthermore, the prediction is that 
for increased neutron number (N = 28) a further rise of the already strong 
deformation of the sulfur isotope chain is to be expected. A similar study 
was performed for heavy nuclei. Fig. 3 presents results for the energy of 
nobelium isotopes, as function of the quadrupole deformation. The hori¬ 
zontal lines mark the experimental binding energy. Measured values for the 
deformation are f32( 252 No) = 0.31F0.02 and /? 2 ( 254 iVo) = 0.32±0.02 u . As 
one can see there is quite a good agreement with experiment. In a survey 
of particle-stable nuclei we performed a scan of nuclear deformation over 
the whole range of possible charge and neutron number. The results are 
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Fig. 2. Quadrupole deformation of neutron-rich sulfur isotopes. Theoretical and exper¬ 
imental 10 values are shown. 
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Fig. 3. Energy of nobelium isotopes as function of quadrupole deformation. The ground 
states of the nuclei exhibit strong prolate deformation.The horizontal lines mark exper¬ 
imental data. 
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shown in Fig. 4. The color-coded information shows the absolute value of 
the quadrupole deformation of the nuclei. Also indicated are the proton 
and neutron drip lines. The results for nuclei beyond the drip line are not 
reliable, as some nucleons start leaking from the nucleus and the numerical 
quantities can depend on the box size of the actual simulation. 

Zooming into the very neutron-rich tin isotopes, one can observe the 
quenching of the shell effects for the Z = 50 nuclei. Fig. 5 shows a pro¬ 
nounced spherical minimum for neutron number N = 82. Increasing the 
number a clear softening of this spherical shape can be observed and slowly 
minima at higher deformation develop. 

As an interesting feature of this plot note the development of a strong 



Fig. 4. Nuclear quadrupole deformation (absolute value) as function of neutron number 
N and charge number Z. The nucleon drip lines are marked. The darker colours mark 
stronger deformation. 


shell closure for the very large neutron number N = 258. This rather 
mo del-independent shell closure was first discussed in 12 . Fig. 6 shows the 
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result for Uranium isotopes for a number of different Skyrme forces and the 
relativistic parameterisation NL-Z2 13 . Here one can observe that already 
for Z = 92 one model shows that the maximum binding in the Uranium 
isotope chain can be obtained for N = 258. Going to higher Z the other 
forces will show a similar behaviour pushing the drip line towards this value. 
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Fig. 5. Binding energy as function of /?2 for various very neutron-rich Sn isotopes. The 
weakening of the effect of the shell closure with increasing neutron number can be clearly 
observed. 

In the usual regime of superheavy nuclei around Z = 114 or Z = 120 
with neutron numbers between N = 170 and N = 190 no clear sign of 
a spherical nucleus can be observed. In addition, the fission barriers are 
quite low of the order of 4 to 5 MeV. In this regime an extended calculation 
allowing for octupole deformations (that have been neglected in the other 
numerical studies shown here) shows a significant reduction of the fission 
barriers as is exemplified in Fig. ??. This general behaviour has also been 
observed in other RMF approaches 14 . However, investigating the corre¬ 
sponding extremely neutron-rich isotopes a different picture emerges. This 
is demonstrated in Fig. 8. This graph shows the deformation-dependent 
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Fig. 6. Binding energies of uranium isotopes for a range of Skyrme forces and a RMF 
parameterisation. 

binding energy of a nucleus with charge Z = 120 and neutron number 
N = 258. In contrast to the solutions of superheavy nuclei with smaller 
neutron numbers and rather small fission barriers, this nucleus exhibits a 
very distinct spherical minimum with a fission barrier of about 9 MeV. A 
similar result was also found in Skyrme models 15 . 

3. Conclusion 

We discussed a general flavor-SU(3) approach to study hadronic and quark 
matter in general, and nuclear structure problems in particular. For a 
quantitatively convincing description of nuclear structure this entails a fine- 
tuning of the parameters 4,6 . 

The resulting model generates good results for binding energies in com¬ 
parison with experimental data. Deformation properties of various nuclei 
also show agreement with experiment. By extending the calculation to very 
neutron-rich nuclei the appearance of a strong shell closure at a (very high) 
neutron number (N = 258) is quite remarkable. This also leads to extremely 
stable superheavy nuclei. Whether such states could be produced, for in- 
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Fig. 7. Example of a superheavy nucleus with Z = 120 and N = 174. The (negative) 
binding energy shows a minimum at large prolate deformation (blue points). However, 
including octupole deformations the fission barrier is reduced significantly (red points). 



Fig. 8. Binding energy as function of deformation for the superheavy nucleus 378 120. 
The curve shows a clear minimum at zero deformation with a large fission barrier. 
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stance in the crust of neutron stars and then ejected by a violent process 

or a star merger, is still completely open and should be estimated. 
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